3.7.99 \(\int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx\) [699]

3.7.99.1 Optimal result
3.7.99.2 Mathematica [B] (warning: unable to verify)
3.7.99.3 Rubi [A] (verified)
3.7.99.4 Maple [F]
3.7.99.5 Fricas [F]
3.7.99.6 Sympy [F]
3.7.99.7 Maxima [F]
3.7.99.8 Giac [F]
3.7.99.9 Mupad [F(-1)]

3.7.99.1 Optimal result

Integrand size = 23, antiderivative size = 313 \[ \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx=-\frac {9 a (a+b \sec (c+d x))^{2/3} \tan (c+d x)}{20 b^2 d}+\frac {3 \sec (c+d x) (a+b \sec (c+d x))^{2/3} \tan (c+d x)}{8 b d}+\frac {\left (18 a^2+25 b^2\right ) \operatorname {AppellF1}\left (\frac {1}{2},\frac {1}{2},-\frac {2}{3},\frac {3}{2},\frac {1}{2} (1-\sec (c+d x)),\frac {b (1-\sec (c+d x))}{a+b}\right ) (a+b \sec (c+d x))^{2/3} \tan (c+d x)}{20 \sqrt {2} b^3 d \sqrt {1+\sec (c+d x)} \left (\frac {a+b \sec (c+d x)}{a+b}\right )^{2/3}}-\frac {a \left (9 a^2+11 b^2\right ) \operatorname {AppellF1}\left (\frac {1}{2},\frac {1}{2},\frac {1}{3},\frac {3}{2},\frac {1}{2} (1-\sec (c+d x)),\frac {b (1-\sec (c+d x))}{a+b}\right ) \sqrt [3]{\frac {a+b \sec (c+d x)}{a+b}} \tan (c+d x)}{10 \sqrt {2} b^3 d \sqrt {1+\sec (c+d x)} \sqrt [3]{a+b \sec (c+d x)}} \]

output
-9/20*a*(a+b*sec(d*x+c))^(2/3)*tan(d*x+c)/b^2/d+3/8*sec(d*x+c)*(a+b*sec(d* 
x+c))^(2/3)*tan(d*x+c)/b/d+1/40*(18*a^2+25*b^2)*AppellF1(1/2,-2/3,1/2,3/2, 
b*(1-sec(d*x+c))/(a+b),1/2-1/2*sec(d*x+c))*(a+b*sec(d*x+c))^(2/3)*tan(d*x+ 
c)/b^3/d/((a+b*sec(d*x+c))/(a+b))^(2/3)*2^(1/2)/(1+sec(d*x+c))^(1/2)-1/20* 
a*(9*a^2+11*b^2)*AppellF1(1/2,1/3,1/2,3/2,b*(1-sec(d*x+c))/(a+b),1/2-1/2*s 
ec(d*x+c))*((a+b*sec(d*x+c))/(a+b))^(1/3)*tan(d*x+c)/b^3/d/(a+b*sec(d*x+c) 
)^(1/3)*2^(1/2)/(1+sec(d*x+c))^(1/2)
 
3.7.99.2 Mathematica [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(19015\) vs. \(2(313)=626\).

Time = 45.82 (sec) , antiderivative size = 19015, normalized size of antiderivative = 60.75 \[ \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx=\text {Result too large to show} \]

input
Integrate[Sec[c + d*x]^4/(a + b*Sec[c + d*x])^(1/3),x]
 
output
Result too large to show
 
3.7.99.3 Rubi [A] (verified)

Time = 1.00 (sec) , antiderivative size = 324, normalized size of antiderivative = 1.04, number of steps used = 13, number of rules used = 12, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.522, Rules used = {3042, 4352, 27, 3042, 4570, 27, 3042, 4495, 3042, 4321, 156, 155}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {\csc \left (c+d x+\frac {\pi }{2}\right )^4}{\sqrt [3]{a+b \csc \left (c+d x+\frac {\pi }{2}\right )}}dx\)

\(\Big \downarrow \) 4352

\(\displaystyle \frac {3 \int \frac {\sec (c+d x) \left (-6 a \sec ^2(c+d x)+5 b \sec (c+d x)+3 a\right )}{3 \sqrt [3]{a+b \sec (c+d x)}}dx}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\int \frac {\sec (c+d x) \left (-6 a \sec ^2(c+d x)+5 b \sec (c+d x)+3 a\right )}{\sqrt [3]{a+b \sec (c+d x)}}dx}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\int \frac {\csc \left (c+d x+\frac {\pi }{2}\right ) \left (-6 a \csc \left (c+d x+\frac {\pi }{2}\right )^2+5 b \csc \left (c+d x+\frac {\pi }{2}\right )+3 a\right )}{\sqrt [3]{a+b \csc \left (c+d x+\frac {\pi }{2}\right )}}dx}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 4570

\(\displaystyle \frac {\frac {3 \int \frac {\sec (c+d x) \left (3 a b+\left (18 a^2+25 b^2\right ) \sec (c+d x)\right )}{3 \sqrt [3]{a+b \sec (c+d x)}}dx}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\frac {\int \frac {\sec (c+d x) \left (3 a b+\left (18 a^2+25 b^2\right ) \sec (c+d x)\right )}{\sqrt [3]{a+b \sec (c+d x)}}dx}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\frac {\int \frac {\csc \left (c+d x+\frac {\pi }{2}\right ) \left (3 a b+\left (18 a^2+25 b^2\right ) \csc \left (c+d x+\frac {\pi }{2}\right )\right )}{\sqrt [3]{a+b \csc \left (c+d x+\frac {\pi }{2}\right )}}dx}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 4495

\(\displaystyle \frac {\frac {\frac {\left (18 a^2+25 b^2\right ) \int \sec (c+d x) (a+b \sec (c+d x))^{2/3}dx}{b}-\frac {2 a \left (9 a^2+11 b^2\right ) \int \frac {\sec (c+d x)}{\sqrt [3]{a+b \sec (c+d x)}}dx}{b}}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\frac {\frac {\left (18 a^2+25 b^2\right ) \int \csc \left (c+d x+\frac {\pi }{2}\right ) \left (a+b \csc \left (c+d x+\frac {\pi }{2}\right )\right )^{2/3}dx}{b}-\frac {2 a \left (9 a^2+11 b^2\right ) \int \frac {\csc \left (c+d x+\frac {\pi }{2}\right )}{\sqrt [3]{a+b \csc \left (c+d x+\frac {\pi }{2}\right )}}dx}{b}}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 4321

\(\displaystyle \frac {\frac {\frac {2 a \left (9 a^2+11 b^2\right ) \tan (c+d x) \int \frac {1}{\sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1} \sqrt [3]{a+b \sec (c+d x)}}d\sec (c+d x)}{b d \sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1}}-\frac {\left (18 a^2+25 b^2\right ) \tan (c+d x) \int \frac {(a+b \sec (c+d x))^{2/3}}{\sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1}}d\sec (c+d x)}{b d \sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1}}}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 156

\(\displaystyle \frac {\frac {\frac {2 a \left (9 a^2+11 b^2\right ) \tan (c+d x) \sqrt [3]{\frac {a+b \sec (c+d x)}{a+b}} \int \frac {1}{\sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1} \sqrt [3]{\frac {a}{a+b}+\frac {b \sec (c+d x)}{a+b}}}d\sec (c+d x)}{b d \sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1} \sqrt [3]{a+b \sec (c+d x)}}-\frac {\left (18 a^2+25 b^2\right ) \tan (c+d x) (a+b \sec (c+d x))^{2/3} \int \frac {\left (\frac {a}{a+b}+\frac {b \sec (c+d x)}{a+b}\right )^{2/3}}{\sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1}}d\sec (c+d x)}{b d \sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)+1} \left (\frac {a+b \sec (c+d x)}{a+b}\right )^{2/3}}}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

\(\Big \downarrow \) 155

\(\displaystyle \frac {\frac {\frac {\sqrt {2} \left (18 a^2+25 b^2\right ) \tan (c+d x) (a+b \sec (c+d x))^{2/3} \operatorname {AppellF1}\left (\frac {1}{2},\frac {1}{2},-\frac {2}{3},\frac {3}{2},\frac {1}{2} (1-\sec (c+d x)),\frac {b (1-\sec (c+d x))}{a+b}\right )}{b d \sqrt {\sec (c+d x)+1} \left (\frac {a+b \sec (c+d x)}{a+b}\right )^{2/3}}-\frac {2 \sqrt {2} a \left (9 a^2+11 b^2\right ) \tan (c+d x) \sqrt [3]{\frac {a+b \sec (c+d x)}{a+b}} \operatorname {AppellF1}\left (\frac {1}{2},\frac {1}{2},\frac {1}{3},\frac {3}{2},\frac {1}{2} (1-\sec (c+d x)),\frac {b (1-\sec (c+d x))}{a+b}\right )}{b d \sqrt {\sec (c+d x)+1} \sqrt [3]{a+b \sec (c+d x)}}}{5 b}-\frac {18 a \tan (c+d x) (a+b \sec (c+d x))^{2/3}}{5 b d}}{8 b}+\frac {3 \tan (c+d x) \sec (c+d x) (a+b \sec (c+d x))^{2/3}}{8 b d}\)

input
Int[Sec[c + d*x]^4/(a + b*Sec[c + d*x])^(1/3),x]
 
output
(3*Sec[c + d*x]*(a + b*Sec[c + d*x])^(2/3)*Tan[c + d*x])/(8*b*d) + ((-18*a 
*(a + b*Sec[c + d*x])^(2/3)*Tan[c + d*x])/(5*b*d) + ((Sqrt[2]*(18*a^2 + 25 
*b^2)*AppellF1[1/2, 1/2, -2/3, 3/2, (1 - Sec[c + d*x])/2, (b*(1 - Sec[c + 
d*x]))/(a + b)]*(a + b*Sec[c + d*x])^(2/3)*Tan[c + d*x])/(b*d*Sqrt[1 + Sec 
[c + d*x]]*((a + b*Sec[c + d*x])/(a + b))^(2/3)) - (2*Sqrt[2]*a*(9*a^2 + 1 
1*b^2)*AppellF1[1/2, 1/2, 1/3, 3/2, (1 - Sec[c + d*x])/2, (b*(1 - Sec[c + 
d*x]))/(a + b)]*((a + b*Sec[c + d*x])/(a + b))^(1/3)*Tan[c + d*x])/(b*d*Sq 
rt[1 + Sec[c + d*x]]*(a + b*Sec[c + d*x])^(1/3)))/(5*b))/(8*b)
 

3.7.99.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 155
Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_)) 
^(p_), x_] :> Simp[((a + b*x)^(m + 1)/(b*(m + 1)*Simplify[b/(b*c - a*d)]^n* 
Simplify[b/(b*e - a*f)]^p))*AppellF1[m + 1, -n, -p, m + 2, (-d)*((a + b*x)/ 
(b*c - a*d)), (-f)*((a + b*x)/(b*e - a*f))], x] /; FreeQ[{a, b, c, d, e, f, 
 m, n, p}, x] &&  !IntegerQ[m] &&  !IntegerQ[n] &&  !IntegerQ[p] && GtQ[Sim 
plify[b/(b*c - a*d)], 0] && GtQ[Simplify[b/(b*e - a*f)], 0] &&  !(GtQ[Simpl 
ify[d/(d*a - c*b)], 0] && GtQ[Simplify[d/(d*e - c*f)], 0] && SimplerQ[c + d 
*x, a + b*x]) &&  !(GtQ[Simplify[f/(f*a - e*b)], 0] && GtQ[Simplify[f/(f*c 
- e*d)], 0] && SimplerQ[e + f*x, a + b*x])
 

rule 156
Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_)) 
^(p_), x_] :> Simp[(e + f*x)^FracPart[p]/(Simplify[b/(b*e - a*f)]^IntPart[p 
]*(b*((e + f*x)/(b*e - a*f)))^FracPart[p])   Int[(a + b*x)^m*(c + d*x)^n*Si 
mp[b*(e/(b*e - a*f)) + b*f*(x/(b*e - a*f)), x]^p, x], x] /; FreeQ[{a, b, c, 
 d, e, f, m, n, p}, x] &&  !IntegerQ[m] &&  !IntegerQ[n] &&  !IntegerQ[p] & 
& GtQ[Simplify[b/(b*c - a*d)], 0] &&  !GtQ[Simplify[b/(b*e - a*f)], 0]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 4321
Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_ 
Symbol] :> Simp[Cot[e + f*x]/(f*Sqrt[1 + Csc[e + f*x]]*Sqrt[1 - Csc[e + f*x 
]])   Subst[Int[(a + b*x)^m/(Sqrt[1 + x]*Sqrt[1 - x]), x], x, Csc[e + f*x]] 
, x] /; FreeQ[{a, b, e, f, m}, x] && NeQ[a^2 - b^2, 0] &&  !IntegerQ[2*m]
 

rule 4352
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + ( 
a_))^(m_), x_Symbol] :> Simp[(-d^3)*Cot[e + f*x]*(a + b*Csc[e + f*x])^(m + 
1)*((d*Csc[e + f*x])^(n - 3)/(b*f*(m + n - 1))), x] + Simp[d^3/(b*(m + n - 
1))   Int[(a + b*Csc[e + f*x])^m*(d*Csc[e + f*x])^(n - 3)*Simp[a*(n - 3) + 
b*(m + n - 2)*Csc[e + f*x] - a*(n - 2)*Csc[e + f*x]^2, x], x], x] /; FreeQ[ 
{a, b, d, e, f, m}, x] && NeQ[a^2 - b^2, 0] && GtQ[n, 3] && (IntegerQ[n] || 
 IntegersQ[2*m, 2*n]) &&  !IGtQ[m, 2]
 

rule 4495
Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(cs 
c[(e_.) + (f_.)*(x_)]*(B_.) + (A_)), x_Symbol] :> Simp[(A*b - a*B)/b   Int[ 
Csc[e + f*x]*(a + b*Csc[e + f*x])^m, x], x] + Simp[B/b   Int[Csc[e + f*x]*( 
a + b*Csc[e + f*x])^(m + 1), x], x] /; FreeQ[{a, b, A, B, e, f, m}, x] && N 
eQ[A*b - a*B, 0] && NeQ[a^2 - b^2, 0]
 

rule 4570
Int[csc[(e_.) + (f_.)*(x_)]*((A_.) + csc[(e_.) + (f_.)*(x_)]*(B_.) + csc[(e 
_.) + (f_.)*(x_)]^2*(C_.))*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_S 
ymbol] :> Simp[(-C)*Cot[e + f*x]*((a + b*Csc[e + f*x])^(m + 1)/(b*f*(m + 2) 
)), x] + Simp[1/(b*(m + 2))   Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m*Simp[ 
b*A*(m + 2) + b*C*(m + 1) + (b*B*(m + 2) - a*C)*Csc[e + f*x], x], x], x] /; 
 FreeQ[{a, b, e, f, A, B, C, m}, x] &&  !LtQ[m, -1]
 
3.7.99.4 Maple [F]

\[\int \frac {\sec \left (d x +c \right )^{4}}{\left (a +b \sec \left (d x +c \right )\right )^{\frac {1}{3}}}d x\]

input
int(sec(d*x+c)^4/(a+b*sec(d*x+c))^(1/3),x)
 
output
int(sec(d*x+c)^4/(a+b*sec(d*x+c))^(1/3),x)
 
3.7.99.5 Fricas [F]

\[ \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx=\int { \frac {\sec \left (d x + c\right )^{4}}{{\left (b \sec \left (d x + c\right ) + a\right )}^{\frac {1}{3}}} \,d x } \]

input
integrate(sec(d*x+c)^4/(a+b*sec(d*x+c))^(1/3),x, algorithm="fricas")
 
output
integral(sec(d*x + c)^4/(b*sec(d*x + c) + a)^(1/3), x)
 
3.7.99.6 Sympy [F]

\[ \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx=\int \frac {\sec ^{4}{\left (c + d x \right )}}{\sqrt [3]{a + b \sec {\left (c + d x \right )}}}\, dx \]

input
integrate(sec(d*x+c)**4/(a+b*sec(d*x+c))**(1/3),x)
 
output
Integral(sec(c + d*x)**4/(a + b*sec(c + d*x))**(1/3), x)
 
3.7.99.7 Maxima [F]

\[ \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx=\int { \frac {\sec \left (d x + c\right )^{4}}{{\left (b \sec \left (d x + c\right ) + a\right )}^{\frac {1}{3}}} \,d x } \]

input
integrate(sec(d*x+c)^4/(a+b*sec(d*x+c))^(1/3),x, algorithm="maxima")
 
output
integrate(sec(d*x + c)^4/(b*sec(d*x + c) + a)^(1/3), x)
 
3.7.99.8 Giac [F]

\[ \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx=\int { \frac {\sec \left (d x + c\right )^{4}}{{\left (b \sec \left (d x + c\right ) + a\right )}^{\frac {1}{3}}} \,d x } \]

input
integrate(sec(d*x+c)^4/(a+b*sec(d*x+c))^(1/3),x, algorithm="giac")
 
output
integrate(sec(d*x + c)^4/(b*sec(d*x + c) + a)^(1/3), x)
 
3.7.99.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\sec ^4(c+d x)}{\sqrt [3]{a+b \sec (c+d x)}} \, dx=\int \frac {1}{{\cos \left (c+d\,x\right )}^4\,{\left (a+\frac {b}{\cos \left (c+d\,x\right )}\right )}^{1/3}} \,d x \]

input
int(1/(cos(c + d*x)^4*(a + b/cos(c + d*x))^(1/3)),x)
 
output
int(1/(cos(c + d*x)^4*(a + b/cos(c + d*x))^(1/3)), x)